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ABSTRACT 

We present high quality optical spectroscopy of the SU UMa-subtype dwarf nova, 
VWHyi taken while the system was in quiescence. An S-wave is executed by the emis- 
sion cores of the Hydrogen Balmer lines and by the emission lines of Hei, Call, Fell 
and He II. Using Doppler tomography we show it originates in the accretion stream-disc 
impact region. The Hen emission is strongly phase-dependent, suggesting it originates 
exclusively within a hot cavity at the initial impact. We map the ionization structure 
of the stream-disc interaction region. One possible interpretation of this is that the 
Balmer hotspot lies downstream of the Hell hotspot in the outer accretion disc, with 
the He I and Can hotspots at intermediate locations between the two. This suggests 
that Balmer emission is suppressed until material has cooled somewhat downstream 
of the impact site and is able to recombine. We favour a phase offset of 0.15 ± 0.04 
between the photometric ephemeris and inferior conjunction of the mass donor. The 
white dwarf contributes significantly to the optical continuum, with broad Balmer 
absorption and narrow MgnA4481 absorption clearly apparent. This latter feature 
yields the gravitational redshift: u grav = 38 ± 21kms _1 , so Mi = 0. 71 toil M o- Tnis 
implies M2 = 0.11 ±0.03Mq and hence the donor is not a brown dwarf. A prominent 
Balmer jump is also observed. We note that the previously accepted system param- 
eters for both VWHyi and WXHyi incoporate an algebraic error, and we provide a 
recalculated M\{q) plane for WXHyi. 
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1 INTRODUCTION 

The close interacting binary, VWHyi, is a cataclysmic vari- 
able (CV), a system in which a Roche-lobe filling low-mass, 
late-type main sequence donor star transfers matter onto a 
white dwarf (WD) primary via an accretion disc. With a 
quiescent visual magnitude of 13.8, VWHyi is the brightest 
example of the SUUMa-type dwarf novae (DNe). These sys- 
tems display two distinct modes of outburst, namely, normal 
DNe outbursts and less frequent larger amplitude super- 
outbursts. These two modes of outburst occur in VWHyi 
typically about every 30 d and 180 d reaching a peak visual 
magnitude of 9.5 and 8.7 respectively. 

The DNe outburst can be understood in the context 
of the disc instability model (DIM) (Osaki 1974). In this 
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model a thermal-viscous limit-cycle instability (Meyer & 
Meyer-Hofmeister 1981), caused by a sensitive temperature 
dependence of opacity in regions of partial ionization in the 
disc, induces a switch between low and high viscosity states, 
and consequently between low and high mass transfer states. 
There remain some questions about the appropriateness of 
this model however (see Lasota 2001, for a review). The su- 
perout burst, distinguishable by a photometric modulation 
known as superhumps, can be explained by a model com- 
bining the thermal instability of the DNe outbursts with a 
tidal instability caused by the action of the donor in a 3:1 
resonance with particle orbits in the accretion disc. This is 
the thermal-tidal instability (TTI) model of Osaki (1989, 
1996). 

As with other SU UMa-type systems, VWHyi has a pe- 
riod that lies below the CV period gap. Its period, P or b = 
1.78 h, was determined photometrically by way of the orbital 
hump observed at quiescence which arises from the phase 
dependent visibility of the accretion stream-disc impact re- 
gion (i.e. the hotspot) throughout the orbit (Vogt 1974; van 
Amerongen et al. 1987). 
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1.1 The hotspot 

Hotspot emission should arise primarily from dissipation of 
kinetic energy as the in-falling stream impacts the disc edge, 
merges with the Keplerian motion of the disc and radiates 
away during the encounter. The details of the stream-disc 
impact region are as yet poorly understood. Study of the 
interaction between the two flows, the stream from L\ and 
the disc itself, may lead to a better understanding of the 
viscosity in accretion discs. 

The flows are highly supersonic, and so the picture of 
the hydrodynamics of the impact point is believed to be 
one which involves the development of shock fronts: one in 
the disc material and one in the incoming stream (Rozyczka 
& Schwarzenberg-Czerny 1987). Emission may represent a 
mixture of accretion stream and disc material compressed 
between these shock surfaces (Marsh et al. 1990) , or emission 
may be seen at two distinct velocities. The smoothed particle 
hydrodynamics (SPH) simulations of Foulkes et al. (2004) 
show two distinct S- waves in simulated trailed spectrograms. 

Other numerical investigations of the stream-disc im- 
pact have been carried out by Armitage & Livio (1996, 
1998), and by Kunze, Speith & Hessman (2001). The lat- 
ter paper shows, through high resolution SPH simulations, 
a substantial vertical deflection of the impacting accretion 
stream after undergoing shock interaction at the hotspot. 
The material is also seen to overflow the disc, re-impacting 
the disc surface at a secondary hotspot close to the circular- 
ization radius at orbital phase 0.5. 

1.2 System parameters 

Despite being an extensively studied system, particularly in 
the UV, the basic parameters of VW Hyi remain uncertain. 
Schoembs & Vogt (1981) put forward an orbital inclination, 
i, of 60° ± 10°, based on the fact that an orbital hump 
is observed in the lightcurve but no eclipse. Schoembs & 
Vogt (1981) also determined component masses by a method 
based upon that of Warner (1973) relating the radial veloc- 
ity semi-amplitude, Kx, of the emission lines and the pro- 
jected rotation velocity of the disc Vd sini to the mass ra- 
tio, q = Ms/Mx. The estimation of K\ using the emission 
lines in CVs has proved, in general, to be unreliable (see e.g. 
Home, Wood & Stiening (1991) for comparison between the 
spectroscopic determination of K\ and the prediction of K\ 
based on the photometric solution in the eclipsing system 
HT Cas). Often the hotspot distorts the apparent radial ve- 
locities as measured from the emission lines, even if only the 
line wings are used, as discussed by Marsh (1998). Schoembs 
& Vogt (1981) report K x = 78 ± 14 km s" 1 from their radial 
velocity curve of Ha. Their tabulated system parameters 
are, however, inconsistent with this value: Mi = 0.63 Mq, 
q = 1/6 and i = 60° give Ki = 57kms _1 . A missing expo- 
nent (4/3) on the (1 + q) term in their equation 7 appears to 
be the source of this discrepancy. It reappears in their sub- 
sequent equations, but has been omitted when calculating 
the orbital inclination in their figures 6 and 7. Based on the 
authors' method we have recalculated Mi(q) for their val- 
ues of Ki . For WX Hyi we show this in Figure [T] Taking the 
authors' assumptions at face value we find Mi = 0-92±g;if 
and q = 0.2lt°°l for VWHyi, and Mi = 1.14±g;|| and 
q = 0.23±o:ol for WXHyi. However, we feel that these val- 
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Figure 1. Mx(q) plane for WXHyi after figure 7 of Schoembs 
& Vogt (1981) for a range of values of their /3 parameter and 
using their value for K\ of 67 ± 6kms _1 . Broken lines indicate 
the range of values of orbital inclination considered by Schoembs 
& Vogt. In order to directly compare with Schoembs & Vogt's 
figure, the abscissa label we use here follows their definition of 
q, i.e. Mx/Ms, this being the reciprocal of the definition we have 
adopted in the rest of this work (M2/M1). 

ues are highly insecure, as the value of /3 particularly may be 
quite different from that assumed, where /3 is approximately 
the ratio of the primary radius to innermost accretion disc 
radius contributing to the line wings. 

For VWHyi, perhaps more credible, by measurement 
of the gravitational redshift from the shift in photospheric 
absorption lines of the white dwarf, Sion et al. (1997) found 
Mi = 0.861q'32 Mq, the weakest part of this work being the 
estimate of the systemic velocity. Warner & Woudt (2002) 
assume that the minimum period observed for dwarf novae 
oscillations (DNOs) in VW Hyi reflects the Kepler period at 
the surface of the white dwarf, Pk- Together with the mass- 
radius relationship for white dwarfs (Nauenberg 1972), they 
calculate a value for Mi of 0.702 Mq. It seems appropriate 
to consider this value as a possible lower limit on Mi as the 
observed minimum DNO period could exceed Pk- Most re- 
cently, from fitting of the J band spectral energy distribution 
of the donor, Mennickent, Diaz & Tappert (2004) suggest it 
to be of spectral type L0, putting M2 ~ 0.06 Mq, and there- 
fore implying a star below the hydrogen-burning minimum 
mass limit, a brown dwarf. If reliable, this would be the first 
detection of the secondary star. 

There exists a relation between the fractional super- 
hump excess, e, defined as (P s h — P or b)/Porb and the mass 
ratio q, where P s h and P or b are the superhump and orbital 
periods respectively. This has been calibrated empirically 
by Patterson et al. (2005). For VWHyi, e = 0.0331 ± 0.0008 
(Patterson 1998), and this indicates q = 0.148 ± 0.004. In 
this work we will assume Mi = 0.711q'26 Mq (see Section 
13.21 , and the value of q derived from the superhump excess, 
q = 0.148 ± 0.004, unless stated otherwise. 

In this paper, we present high quality optical spec- 
troscopy of VWHyi. The observations and data reduction 
are described in Section 2. In Section 3 we present our re- 
sults, including Doppler maps of VWHyi and determination 
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Table 1. 2003 February 14, CTIO observations of VWHyi 



HJD-2452680 


Exposure 


Exposure 


N 


Orbital 


(mid-exposure) 




Time/s 




Phase 


4.5383 


EG21 


5.0 


1 




4.5404 - 4.5416 


EG21 


20.0 


3 




4.5589 - 4.5831 


VWHyi 


120.0 


14 


0.677-0.003 


4.5887 - 4.6220 


VWHyi 


120.0 


20 


0.079-0.528 



of the gravitational redshift of the white dwarf. Interpreta- 
tion of these results is given in Section 4, which includes 
discussion of the ionization structure of the stream-disc im- 
pact and further discussion of the system parameters for 
VWHyi. A summary is given in Section 5. 



2 OBSERVATIONS AND DATA REDUCTION 

We observed VWHyi on 2003 February 14 01:09-02:56 UT 
using the CTIO 4.0 m R-C spectrograph as detailed in Ta- 
ble The AAVSO lightcurve for VWHyi shows an outburst 
some 27 days before our observations, peaking on 2003 Jan- 
uary 18. The timing of the subsequent outburst is unclear. 
An uncertain magnitude of 9.5 was recorded on 2003 Febru- 
ary 16, while the following day the magnitude was confirmed 
to be 10.6, roughly three magnitudes above quiescence. This 
suggests that our observations took place just prior to the 
onset of outburst. The previous superout burst took place in 
November 2002 (Waagen 2004). 

During the observations the seeing was 1.4". Using 
the KPGL1 grating and a Loral 3K CCD in conjunction 
with a 1.3" slit, a wavelength coverage of 3590-6686 A and 
a FWHM resolution of ~ 3 A as measured from arc lines 
was obtained. We obtained 34 useful spectra, each with 
an integration time of 120 s, together with three approx- 
imately equally interspersed HeNeAr arc-lamp exposures, 
giving ~ 85 per cent orbital coverage. The spectra have a 
signal-to- noise ratio of ~ 50. For the purposes of flux cali- 
bration, the standard star EG21 was observed. 

The observations were reduced following standard pro- 
cedures using IRAF \ The frames were first bias-subtracted, 
flat-fielded and sky-subtracted. The spectra were then ex- 
tracted according to the optimal algorithm of Home (1986). 
The arcs were extracted at the position of the corresponding 
object spectrum which were then used to wavelength cali- 
brate the data by interpolation between neighbouring arcs. 
A zero-point correction of +4.5 ± 1.7 km s was applied to 
the wavelength scale, derived from the O I A5577 sky emis- 
sion line. The spectra were subsequently flux calibrated. All 
spectra were transformed to the heliocentric restframe. 

Without correcting for slit losses, we convolved the flux 
calibrated average spectrum with the filter bandpass, hence 
deducing V < 14.5. This indicates the system was quiescent. 
The system may have been fainter than the usual quiescent 
magnitude, but slit losses of 0.7 magnitudes are possible. 

1 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



3 RESULTS 

The average of all 34 spectra is presented in Figure H and 
is quite remarkable in its detail and wealth of features. The 
spectrum shows strong Balmer lines and a prominent Balmer 
jump. He I, Fell and Call emission features are also present, 
and He II emission is weak but clearly detected. Ha is domi- 
nated by double-peaked emission: a typical signature of disc 
emission, due to the Doppler shifting of orbiting disc mate- 
rial. 

The emission cores of the other lines in the Balmer se- 
ries are flanked by broad absorption wings. We attribute 
this broad absorption to Stark broadening in the photo- 
sphere of the white dwarf, as the line width of this absorption 
would seem to preclude the possibility of kinematic broad- 
ening in the disc. A conservative estimate for the base-to- 
base width of the H7 absorption is 210 A. If this is due to 
Doppler broadening, a radial velocity v T — 7300 kms -1 is re- 
quired. Assuming i — 60°, the corresponding true velocity is 
i>truc — 8400 kms -1 . At the inner disc radius, the Keplerian 
velocity of disc material is at most ~ 4200 kms -1 , assuming 
the highest estimate for Mi , the inner disc radius to coincide 
with the white dwarf surface, and the mass-radius relation 
for white dwarfs of Nauenberg (1972). In CVs, generally 
the accretion flow dominates at optical wavelengths, but it 
is not too unusual to see optical absorption features from 
the white dwarf as we do here, in systems where the mass- 
accretion rate is low for some reason (e.g. Marsh, Home & 
Shipman 1987; Hessman et al. 1989). These quiescent broad 
absorption features are seen in one previous spectral study 
of VWHyi (Hassall et al. 1983), but are in contrast to the re- 
sults of other studies also conducted in quiescence, where no 
such features are observed (Schoembs & Vogt 1981; Saygac 
1989). This might suggest the system being in an unusu- 
ally low state. However there is no clear indication of this 
from the AAVSO lightcurve. In VWHyi, the white dwarf is 
clearly revealed in quiescence in the UV (Mateo & Szkody 
1984; Sion et al. 1995). The Balmer jump shows absorption 
in the Balmer continuum (i.e. A < 3646 A). We refer to 
Gansicke, Beuermann & Thomas (1997) 's spectrum of EK 
TrA, which shows a similar feature although less pronounced 
than that seen in our data. 

The double-peaked emission component, clear in Ha, 
becomes progressively weaker towards the higher order 
Balmer lines, as seen in the trailed spectrograms in the top 
panels of Figure [3] Here, the spectra have been normalized 
using a low-order spline fit to the continuum. We also see 
the S-waves executed by the emission cores, due to the pe- 
riodic redshifting and blueshifting of an isolated emission 
component orbiting in the binary system. The amplitude 
of the S-wave is ~ 500 kms -1 . The emission appears to be 
strongest in the transition from red to blue. Away from the 
main S-wave, the double-peaked component of the emission 
is still variable. The enhanced regions of this double-peaked 
emission may be independent, but their phasing is roughly 
consistent with the sinusoidal motion of a single emission 
site in the disc, as indicated by the second overplotted sine 
wave on the trailed spectrum of Ha. 

S-waves are clearly detected in the Balmer series almost 
right up to the Balmer limit (Smith 2006). The majority of 
the emission lines in the other species also show evidence of 
S-wave emission (see Figure [5] and Sections 13. 51 and 13. 61 . 
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Figure 2. Average spectrum of VWHyi, in which each individual spectrum is given equal weighting and no correction for orbital motion 
has been made. Overplotted is a model white dwarf spectrum (solid blue line) (see Section 13.31 . 



At the suggestion of the referee, and for the purposes 
of comparison with previous studies, we employed the dou- 
ble Gaussian convolution method (Schneider & Young 1980; 
Shaffer 1983) to examine the behaviour of the Ha emission 
profile. Results of this are to be included in Smith (2006). 
The value of K em so obtained is in agreement with that 
expected from our chosen system parameters, but there is 
a strong likelihood of hotspot contamination even if large 
Gaussian separations are used, particularly here where the 
S-wave clearly dominates the disc emission (Figure 0. 



3.1 Doppler tomography 

The indirect imaging technique of Doppler tomography 
(Marsh & Home 1988) translates orbital phase-resolved 
emission line profiles into a distribution of line emission in 
two-dimensional velocity space. In this work, Doppler maps 
were produced using Tom Marsh's MOLLY and DOPPLER 
software. With the exception of Figure|S] the Fourier-filtered 
back-projection (FFBP) method was used (see Marsh 2001, 
for a recent review) . As the exposure time used for all spec- 
tra is 120 s, we can expect a point emission source to be 
blurred by 6.7° azimuthally, which equates to a phase reso- 



lution of 0.019. Radial blurring will be the result of the finite 
spectral resolution of the data (~ 170 kms -1 ). 

Modulation is apparent in the continuum level of the 
spectra, but it is not possible to distinguish real intrinsic 
modulation from that which can be attributed to variable 
slit losses. Hence, the examination of orbital variation in line 
flux is inappropriate. Instead orbital variations of equivalent 
width in Ha and H/3 are plotted (Figure 0J. The underlying 
white dwarf absorption makes the calculation of true equiv- 
alent width somewhat problematic, and we did not attempt 
this for the higher order Balmer lines. Nonetheless the rel- 
ative values for Ha and H/3 should be reliable. The figure 
suggests a weakening of Ha from approximately phase 0.68 
to phase 0.91. Without coverage over multiple orbits, it is 
uncertain whether this is a systematic variation or due to 
flickering. One of the underpinning assumptions of Doppler 
tomography is that line flux is constant over the orbit. To 
fulfill this we normalized the spectra for each map to con- 
stant integrated line flux, except where the effects of noise 
then produced significant artefacts in the maps, in H<5, Hell, 
Can and Fen. 

The S-wave in the trailed spectrograms maps into a con- 
centrated region of emission in velocity space. We estimated 
the systemic velocity for the Doppler maps in two ways. The 



© 2006 RAS, MNRAS OOO.ITHTuI 



VWHyi: optical spectroscopy and Doppler tomography 5 




- l 500- * 000 - 500 500 I0O0 * IO0 500- '000-500 500 '000 ' 500- * 500- 1 000 - 500 500 1000 1500- '500- 1000-SOG 500 'GOO '500 



V. (km s"') V, {um s"') V, (km s*') V, (km s" 1 ) 

Figure 3. Top: Trailed spectra of Ha, H/3, H7 and H<5 emission respectively, where the spectra have been averaged into 20 orbital phase 
bins and the full cycle repeated for clarity. Overplotted on Ha are sine waves tracing the motion of the main S-wave emission (black), 
and indicating the motion of a second possible isolated emission site (dark orange). Bottom: Doppler maps of Balmer emission. The 
system parameters used in the construction of the maps are shown in the legends. The value for i given by Schoembs & Vogt (1981) was 
assumed. <f>o gives the orbital phase with respect to the photometric ephemeris of van Amerongen et al. (1987). The ballistic trajectory of 
the gas stream is represented by the solid line, with open circles denoting steps of 0.2Rli along the stream. The 'Kepler shadow', which 
shows the Keplerian velocity at equivalent positions along the ballistic trajectory, is represented by the dotted line with filled circles. 
The crosses, reading from top to bottom, denote the centre of mass of the secondary, the centre of mass of the system, and the centre 
of mass of the primary respectively. The Roche lobe of the primary is given by the dashed line, and the Roche lobe of the secondary is 
shown by a solid line. The circular solid line represents Keplerian velocity at the tidal truncation radius of the accretion disc. 
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Figure 4. Orbital variation of line equivalent width for Ha and 
H/3. 



first was to make the sharpest possible maps, i.e. to use the 
value of 7 which minimised the width of a Gaussian fitted 
to the concentrated emission spot in the Doppler map. The 



average 7 value found by this method for the Balmer lines 
and HelA5876 is 29kms~ 1 , whilst the spread in values was 
18 to 46kms _1 . Fitting to the entire map as opposed to 
isolating just the emission spot yielded an average value of 
56kms _1 , whilst a search for the 7 value which maximized 
the strength of the emission gave 41 km s~ . We adopted 
7 = 40 ± 20kms _1 for the maps. The maps produced using 
any of these differing 7 values are essentially indistinguish- 
able. 



3.2 Gravitational redshift of the white dwarf 

Sion et al. (1997) reported gravitational redshift of 58 ± 
33kms _1 after estimating 7 = 2 ± 14kms _1 , from the ra- 
dial velocity curve of Schoembs & Vogt (1981); Schoembs 
& Vogt (1981) however state that their 7 value is arbi- 
trary and chosen so that the mean radial velocity is zero. 
The narrow absorption feature in the region of He I A4472 is 
MgllA4481 and originates in the white dwarf. In the aver- 
age spectrum this feature has a shift of 78 ± 8kms _1 . We 
can use the difference between the emission line systemic 
velocity, 40 ± 20kms _1 , and the white dwarf MgllA4481 
absorption line, to estimate the white dwarf's gravitational 



redshift. We obtain v s 



38 ±21 kms . This is consistent 



with Sion et al. (1997)'s result and is more secure since it 
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is derived relative to simultaneously measured Balmer emis- 
sion from the low-gravity environment of the outer disc. Us- 
ing this value of fgrav , together with the mass-radius rela- 
tion for a 20 000 K white dwarf with a carbon core (Wood 
1995), we find Mi = 0.71±o.26 M , a white dwarf radius, 
Ri = 9.3^2 e x iO 8 cm, and surface gravity, logg = 8.04±°;|f . 

3.3 Source of the S-wave 

van Amerongen et al. (1987) provide a photometric 
ephemeris of maximum light from the hotspot; the phasing 
of conjunction is unknown. Hence there are two possible sce- 
narios for the origin of the S-wave: the stream-disc impact re- 
gion; or the irradiated face of the donor. Neither hypothesis 
can be immediately ruled out. By varying the system param- 
eters for the map of H7, the likelihood of each was investi- 
gated. Only in the most extreme case that Mwd = 1.44 Mq, 
i.e. the Chandrasekhar mass, could emission correspond to 
the inner face of the donor. Sion (1999) finds the mean mass 
of 76 CV white dwarfs to be 0.86 M , while Webbink (1990) 
finds a geometric mean primary mass for systems below the 
period gap to be 0.61 Mq. A single object of course does not 
have to adhere to such averages, but ~ 1.4 Mq would be 
atypically large. Furthermore, a larger inclination (2 = 75°) 
than estimated for VWHyi and an ultra-low mass ratio 
(0.04) would be required. On the other hand, as Figure El 
shows, the S-wave emission can be naturally attributed to 
the stream-disc impact region using parameters consistent 
with other studies of VWHyi. The coincidence, to within 
0.01 in phase, of the emission spot with the photometric 
hotspot ephemeris (i.e. when no phase offset is applied) is 
also persuasive. 

It is not possible to locate the position of the hotspot 
exactly. It would be expected to lie somewhere in the re- 
gion between the ballistic stream trajectory and its Ke- 
pler shadow. A mixing of the stream and disc material and 
their respective velocities is expected as the material is com- 
pressed between two shock surfaces where the in-falling ac- 
cretion stream impacts the disc edge and merges with the 
Keplerian flow (Marsh et al. 1990). An orbital phase offset 
of 0.15 ±0.04 (where the error comes from extrapolating the 
epochs) with respect to the photometric ephemeris of van 
Amerongen et al. (1987) was required to place the hotspot 
in H7 on the ballistic stream trajectory close to the intersec- 
tion with the tidal radius. This phase offset is comparable 
with that between <j>o and the hotspot orbital hump for the 
similar, but eclipsing, system, ZCha, of 0.18 (Warner 1974). 

3.4 Balmer emission 

Figure |3 shows the Doppler maps in the Balmer lines, con- 
structed using the primary mass found in Section [3.21 com- 
bined with q calculated from the superhump excess, and us- 
ing the phase offset found above. In order to suppress high- 
frequency noise, a Gaussian cutoff is included in the Fourier 
filter that is applied to the spectra prior to back-projection. 
This Gaussian is chosen so that a suitable trade-off between 
resolution and noise is obtained. Here, the FWHM of the 
Gaussian window function is 10, and so there is minimal 
reduction in noise to preserve resolution. 

There is a potential problem with contamination by 



the white dwarf absorption features, as they would be ex- 
pected to have a different (unknown) phasing to that of the 
(hotspot) emission. To account for this, a model white dwarf 
spectrum was constructed, based on the properties found for 
the white dwarf in VWHyi by Godon et al. (2004) from UV 
observations. This model was scaled in flux to match our 
data at the blue end (Figure and then fitted to the ab- 
sorption wings of the Balmer line profiles in the individual 
spectra by allowing shifts in wavelength. The fitted white 
dwarf model spectrum was then subtracted. H7 has the 
cleanest absorption wings in our data and so these were 
used in the fitting, although using the wings of H5 produced 
very similar results. The wavelength shifts obtained in this 
way were not inconsistent with the predicted motion of the 
white dwarf, if the parameters we have used in producing 
the Doppler maps are assumed (Smith 2006). It was found 
that most of the structure seen in the Doppler maps was 
unaffected by this procedure, the hotspot emission particu- 
larly so. What did change was the appearance of a confined 
central emission feature in the region of the white dwarf. 
The central regions of the lines corresponding to this low 
velocity emission will be where the effect of the correction 
is greatest owing to the deep absorption cores of the white 
dwarf spectrum. This feature must, therefore, be considered 
an artefact due to inaccuracies in the fitting procedure and 
uncertainties in the appropriateness of the parameters used 
in constructing the model spectrum. We nevertheless con- 
clude that the effect on the Doppler maps of the white dwarf 
absorption is very minimal and does not need to be consid- 
ered further. 

The maps reveal how the systematic variation through 
lines of the Balmer series seen in the trailed spectrograms 
relates to distribution in velocity space, from a smoother 
distribution over the disc in Ha, to a concentrated hotspot 
in US. 

The overplotted sine waves in the trailed spectrum of 
Ha in Figure [3] indicated that the variable double-peaked 
emission may be correlated. A possible source for this emis- 
sion is a secondary region of enhanced emission seen in Fig- 
ure Here the maximum entropy method (MEM) is used 
(discussed in detail in Marsh & Home (1988)) as opposed 
to FFBP. Both methods produce qualitatively very similar 
maps, depending somewhat on the value of the controlling 
parameter used (FWHM of noise filter in the case of FFBP, 
and x 2 m the case of MEM), but MEM appears to better 
detect detail such as a secondary emission region as is the 
case here. 

3.5 He 1 and He 11 emission 

Figure |S| shows the trailed spectra and corresponding 
Doppler maps for He I A5876 and He I A4472. The Mgll A4481 
absorption feature from the white dwarf creates a dark ring 
artefact in the He I A4472 map. In both lines there is little or 
no evidence for a disc component in emission, while a strong 
S-wave is exhibited, again originating in the stream-disc im- 
pact region. The He I A5016 and He I A4922 lines are possibly 
contaminated by blending with FellA5018 and FellA4924 
respectively, and Doppler maps for these lines show a very 
much less distinct hotspot. Of particular note is the presence 
of He 11 A4686 in the spectra. This also demonstrates S-wave 
motion, with the emission, in contrast to that of the Balmer 
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Figure 5. Doppler map of Ha emission using the maximum en- 
tropy method. Symbols are as in Figure 151 



lines, strongest in the transition from blue to red. This S- 
wave also maps to approximately this same hotspot region 
(Figure |SJ indicating a high temperature at the hotspot. 
Because HellA4686 is a significantly weaker feature in the 
spectra, noise posed a greater problem in the production of 
the Doppler map. Here a Gaussian window function with a 
FWHM of 0.3 was used with the inevitable loss of resolution, 
but still the hotspot remains clear. 

That Hen is detected and seen to emanate from the 
hotspot region in a quiescent DN is relatively rare. It was 
seen in U Gem (Marsh et al. 1990; Groot 2001), where the 
velocity of emitting material was found to lie intermedi- 
ate between that of the accretion stream and the Keplerian 
disc at the stream-impact, and in the newly identified cata- 
clysmic variable Ha0242-2802 (Mason & Howell 2005). 



3.6 Emission in metal lines 

Figure Q shows the corresponding trailed spectra and 
Doppler maps for CallK and FellA5169 respectively. 
FellA5169 may be blended with Mgl, blurring the hotspot. 
The CallK hotspot is clear. There is a suggestion of weak 
disc emission in both Can and Fen. 
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Figure 7. Top left and top right panels show trailed spectra of 
Ca II K and Fe II A5169 respectively, and below the corresponding 
Doppler maps. Symbols are as in Figure O 



maps of other quiescent dwarf novae (Skidmore et al. 2000; 
Rolfe, Abbott & Haswell 2001), and in the quiescent soft X- 
ray transient A0620-00 (Shahbaz et al. 2004). In the latter 
system, a significant modulation of Ha with orbital phase 
is observed, perhaps as a result of a changing perspective of 
an optically thick, elongated hotspot. 

In Figure a second region of enhanced emission in Ha 
is seen. The limited phase coverage of our observations lim- 
its confidence in the reality of this feature. The dark band 
across the maximum-entropy map running from bottom left 
to top right for example, corresponds to the phase gap in 
the observations, which could also influence the appearance 
of this second enhanced emission feature. We note that the 
phase at which the presence of a putative secondary S-wave 
is most convincing corresponds to when Ha emission was 
strongest (Figure |1J. If real, then it would indeed be in- 
teresting, and we are without an explanation for its pres- 
ence. It does not fall where an overflowing accretion stream 
would be expected to re-impact the accretion disc, cf. Kunze 
et al. (2001) and their figure 13. Livio (1993, and references 
therein) discuss secondary interactions between the accre- 
tion stream and disc, and observation of multiple enhanced 
regions of emission can be found in Marsh & Home (1990); 
Neustroev (1998, 2002); Roelofs et al. (2005). 



4 INTERPRETATION 

4.1 The Balmer Doppler maps 

The Doppler maps of the Balmer lines show a progression 
from hotspot emission superimposed on disc emission in 
Ha to a relatively more dominant hotspot in H/3, to solely 
hotspot emission in H5. This may be interpreted as the 
hotspot being optically thick in Ha and less so in H/3, so 
that in Ha the hotspot is weaker relative to the optically 
thin disc emission. This behaviour is seen in the Doppler 



4.2 Ionization structure of the accretion 
stream-disc impact region 

Pringle et al. (1987) estimated the colour temperature of 
the orbital hump due to the hotspot in VWHyi to be 
11 000 ± 2000 K, which is typical of the blackbody temper- 
ature of the hotspot in DNe. Our observations show that 
temperatures in the stream-disc interaction region must be 
sufficient to ionize helium. We show the locations of the 
centroid of concentrated hotspot emission (with <fio = 0.15 
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Figure 6. Top panels left to right show the trailed spectra of HelA4472, HelA5876 and Hen A4686 respectively, and below the corre- 
sponding Doppler maps. Symbols are as in Figure 131 



assumed) in Figure |5] As the true phasing of the system is 
unknown, these positions may be rotated around the centre 
of mass of the system, but their relative positions will re- 
main as shown. The hotspots in the Balmer and He I lines 
lie downstream of the hotspot in HellA4686. This situation 
would be expected if the temperature of the hotspot region 
is such that neutral emission is suppressed by ionization at 
the impact, and that recombination can only take place once 
the stream material mixing with the disc material has moved 
downstream of the initial shock interaction and cooled suf- 
ficiently. As the first ionization potential for Ca is less than 
that of He, so we can also expect Ca II K hotspot emission to 
be seen a lower temperatures downstream of the He II emis- 
sion, just as Figure |H| indicates. From eclipse measurements 
in WZSge, Spruit & Rutten (1998) found the hotspot in 
the optical continuum and in Ha to be separated in phase 
by 0.015, with the continuum leading Ha. In IYUMa, emis- 
sion in the Balmer lines is found to be weaker in the disc 
regions around the stream-disc impact (Rolfe et al. 2005), 
indicating that Balmer emission is indeed suppressed at the 
high temperatures of the hotspot. We similarly see a sup- 
pression of disc Ha emission along the Kepler shadow of the 
stream (Figures [3] and Q)J. Spruit & Rutten (1998) invoke a 
small, fully ionized 'Stromgren sphere' around the hotspot 
to explain their observations of WZ Sge: intense UV emis- 
sion causing a region of reduced Ha brightness in the disc 
bordering the impact point. 

However, the same physical position can map to differ- 
ing velocity locations, and it is also possible that the hotspot 
emission in the different lines emanates from the same phys- 
ical point but that the dynamics of emitting material dif- 
fers. Another possible interpretation of Figure |H| therefore, 



is reached by noting that the Hell hotspot is closer to the 
Kepler shadow than the Can, He I and Balmer hotspots. 
This might indicate that He II emission may arise from ma- 
terial with predominately disc-like velocities, the emission 
perhaps taking place after a greater degree mixing with the 
disc material, while the Balmer emission is dominated by 
material in the ballistic stream. Note 4>o = 0.15 was found 
by placing the H7 hotspot on the stream trajectory. If we 
had instead placed the He 11 hotspot on the stream trajectory 
(effectively rotating the points in Figure |H] about the centre 
of mass, while leaving the stream and Roche lobe loci fixed) 
then the hotspot emission in the other species would fall be- 
low the stream trajectory in the map. In eclipsing systems 
where the phasing is known (e.g. IYUMa, Rolfe et al. 2005) 
the Balmer and He I hotspots occur between the stream tra- 
jectory and its Kepler shadow as they do with our choice of 

00- 

The differing velocities of the hotspot in the various 
spectral lines shows this localized emission does emanate 
from the stream-disc impact. An origin on the mass donor 
star would require all species to have velocities consistent 
with the donor star. Figure |S| shows, however, the spread 
in velocities is bigger than the velocity image of the donor 
Roche lobe. 

Emission in Hell, in contrast to He I, shows a strong 
phase dependence (Figure |SJ (the apparent phase depen- 
dence in HelA4472 we attribute to the MgllA4481 absorp- 
tion). Hell emission requires high temperatures which, in 
quiescence, will likely only be reached at the initial shock- 
heated stream-disc impact point, hence Hell emission is ex- 
pected to be highly localized. Steeghs et al. (2003) discovered 
an unusual feature in their orbital lightcurve of IY UMa: a 
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Figure 8. Detail of hotspot region in velocity space with the 
centroid locations of the hotspots in the emission lines shown. 
Symbols used are indicated in the legend. Additional symbols 
used are as in Figure |3 The dashed line connects the ballistic 
stream and Keplerian disc velocities at a radius of 0.5 f° r the 
system parameters assumed. 



rise occurring between the end of white dwarf ingress and 
hotpot ingress. This is best explained by the brief view into a 
shock-heated cavity with a narrow opening angle, caused by 
the impact of the accretion stream. A similar cavity, though 
wider (as VWHyi has a lower inclination), could cause the 
phase dependence of Hell. The emission might be expected 
to be seen at phases where the hotspot faces us, reaching a 
maximum at the point where we are looking into the cavity 
near maximum redshift, as our data shows. He I emission, on 
the other hand, probably comes from a larger physical vol- 
ume at lower temperatures, and is hence visible throughout 
the orbit. 



4.3 System parameters revisited 

The radius of the accretion disc is unknown. Because of this 
uncertainty it was not possible to discriminate decisively 
between the various system parameters put forward in the 
literature, and all, including the Mi value we derive in Sec- 
tion OH appear to be consistent with the hypothesis of the 
hotspot as being the source of the S-wave observed (Figure 
Ej. From Figure 00 it might appear that the hotspot in sev- 
eral lines falls outside the tidal radius (the disc is inverted in 
velocity space) . Under the assumption of Keplerian disc ve- 
locities, the requirement for the position of the H7 hotspot 
to be within the tidal radius places some tentative limits on 
the system parameters for VWHyi. If q — 0.148 is assumed, 
then Mi sin 3 i < 0.39 ± 0.05 Mq is required. However, we 
may expect the velocities in the outer disc to be ~ 10 % 
sub-Keplerian for a system of this mass-ratio (Truss, pri- 
vate communication). This, coupled with the blurring in 
the Doppler maps, removes the apparent conflict between 
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Figure 9. Doppler maps of H7 showing the location of the 
hotspot using various parameters in the literature, and assuming 
(f>0 = 0.15, as shown in the legend of each panel. Again, symbols 
are as in Figure 13 The Doppler map in the left panel uses the 
value for Mi given by Sion et al. (1997), where q = 0.148 is as- 
sumed, while that in the right panel uses the values of Mennickent 
ct al. (2004). 



the hotspot positions and tidal radius in Figure 03 We note, 
though, that the disc is expected to reside well within -Rtidos 
when in a quiescent state. 



5 SUMMARY 

The following summarises our main findings: 

• We observed VWHyi in quiescence. Slit losses mean 
V < 14.5. 

• There is strong S-wave emission in the cores of the 
Balmer series and in He I, Hell and Call. This S-wave arises 
in the stream-disc impact region. 

• We favour a phase offset of <j)o = 0.15 ± 0.04 between 
the photometric ephemeris and inferior conjunction of the 
mass donor. 

• A progression in the appearance of the Doppler maps, 
from hotspot emission superimposed on disc emission in Ha, 
to relatively more prominent hotspot emission in the higher 
order Balmer lines, is seen. This may be interpreted as the 
hotspot being optically thick in Ha and less so in H/3, so 
that in Ha the hotspot is weaker relative to the optically 
thin disc emission. 

• The stream-disc interaction region shows an ioniza- 
tion structure in velocity space. One interpretation is that 
the Balmer and He I hotspots lie downstream of the Hell 
hotspot. This may be expected if immediate post-shock tem- 
peratures are such that H and He are ionized and recombi- 
nation can only occur once material has cooled downstream 
of the impact point. 

• He 11 emission shows a strong phase dependence, which 
may be the result of restricted viewing of a localized emission 
site in a shock-heated cavity at the stream-disc impact site. 

• A secondary region of enhanced emission is detected in 
the Doppler map of Ha produced by the maximum entropy 
method, though it is unclear how this might arise. 

• We detect broad Balmer absorption and narrow 
MgnA4481 absorption from the white dwarf. We use this 
MgnA4481 feature to determine the gravitational redshift, 
!)„ av = 38 ± 21 kms -1 . 
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• Our preferred system parameters are M 1 = 
(mioieMo (from our u grav ) and q = 0.148 ± 0.004 
from the superhump period excess (Patterson 1998). This 
implies M2 = 0.11 ± 0.03 Mq. Hence the donor star would 
not be a brown dwarf. 

Though Hen in the hotspot of quiescent DN has been 
rarely detected before now, its presence may actually be a 
widespread phenomenon. The S/N of existing observations 
may simply not be sufficiently high to detect it. Our results 
show that modest amounts of 4 m telescope time (less than 
2 hours in this case) can resolve uncertainties regarding the 
fundamental parameters of CVs. 
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